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Abstract 

The  performance  of  a  molten  carbonate  fuel  cell  stack  with  regard  to  safe  and  efficient  electricity  generation  has  been 
investigated  using  the  computational  fluid  dynamics  (CFD)  technique.  A  numerical  model  is  developed,  and  it  is  employed 
to  calculate  the  three-dimensional  distributions  of  the  crucial  parameters  (e.g.,  temperature,  pressure,  concentration,  and 
density)  across  a  stack.  In  particular,  the  model  can  consider  simultaneously  the  dominant  processes  of  a  stack,  such  as  mass 
transport,  chemical  reactions,  heat  transfer,  and  the  voltage-current  relations.  Moreover,  it  is  also  capable  of  calculating  the 
mass  distribution  across  the  stack  rather  than  assuming  a  uniform  distribution.  In  this  paper,  the  model  is  demonstrated  by 
applying  it  to  calculate  fuel  cells  with  three  different  manifolds  (e.g.,  co-,  counter-  and  cross-flow)  in  a  stack.  The  model  is 
an  effective  numerical  tool  for  optimal  design  and  operational  analysis  of  fuel-cell  stacks,  e.g.  for  comparing  performance  with 
different  manifolds  and  to  identify  operational  characteristics. 
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1.  Introduction 

Most  experimental  investigations  of  fuel  cells  are 
expensive.  In  addition,  they  can  only  be  carried  out 
on  a  small  scale,  with  a  few  designs  and  with  difficulties 
of  measurement,  e.g.,  of  the  electrical  resistance.  On 
the  other  hand,  the  cost  of  numerical  modelling  of  fuel 
cells  is  low,  and  it  is  easy  to  set  up  a  number  of 
alternative  new  designs.  With  a  model,  the  quantitative 
prediction  of  the  overall  performance  of  fuel  cells  can 
be  made  over  a  wide  range  of  operations  and  for 
different  transient  conditions.  Hence,  the  numerical 
simulation  of  fuel  cells  is  worthwhile  and  has  been 
used  extensively  [1]. 

The  numerical  simulation  can  be  appropriately  di¬ 
vided  into  modelling  at  the  electrode  level  (micro¬ 
modelling),  at  cell  level,  and  at  stack  level  (macro¬ 
modelling).  Micromodelling  aims  at  building  better 
electrodes  through  study  of  microscopic  processes.  How¬ 
ever,  macromodelling  aims  at  optimizing  design  alter¬ 
natives  and  determining  operational  strategies,  provided 
that  the  electrochemical  performance  of  the  given  elec¬ 
trodes  and  electrolytes  are  known.  In  this  paper,  ma¬ 
cromodelling  of  fuel  cells  is  emphasized.  In  addition, 
Dutch  fuel-cell  research  and  demonstration  programs 
have  a  special  interest  in  the  molten  carbonate  fuel 
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cells  (MCFC)  [2].  Thus,  this  paper  is  exclusively  aimed 
at  MCFC  modelling.  Without  further  specification,  the 
term  ‘fuel  cell’  refers  to  ‘MCFC’  throughout  the  paper. 

In  macromodelling  of  the  MCFC,  the  two-dimensional 
numerical  model  developed  by  Wolf  and  Wilemski  [3] 
is  one  of  the  most  sophisticated  fuel-cell  models  at  the 
cell  level.  However,  there  are  still  few  models  at  the 
stack  level,  and  an  equivalent  cell  is  normally  used  to 
represent  stack  performance.  In  the  parallel  field  of 
macromodelling  of  the  solid  oxide  fuel  cell  (SOFC), 
the  latest  stack  model  introduced  by  Achenbach  [4] 
has  considered  the  three-dimensional  current  and  tem¬ 
perature  distribution  across  the  stack,  but  it  is  still 
based  on  the  assumption  of  a  uniform  mass  distribution 
along  the  stack  direction. 

On  the  other  hand,  experimental  results  recently 
reported  by  Takashima  [5]  have  indicated  that  there 
is  a  significant  inhomogeneity  (e.g.  mass  flow  rate)  along 
the  direction  of  the  stack,  and  it  increases  with  stack 
height  and  operating  pressure.  This  indicates  that  using 
a  two-dimensional  cell  model  to  represent  the  stack 
or  assuming  uniform  mass  distribution  along  the  stack 
direction  is  no  longer  appropriate  for  high  stacks  or 
for  stacks  operating  at  high  pressure. 

Therefore,  an  effective  stack  model  for  MCFC  should 
be  able  to  deal  with  the  physical  and  chemical  phe- 
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nomena  across  the  stack,  and  capable  of  calculating 
the  mass  distribution  along  the  stack  direction.  Cor¬ 
responding  to  these  requirements,  development  of  a 
three-dimensional  simulation  model  is  required  to  take 
account  for  these  effects.  In  this  paper,  the  compu¬ 
tational  fluid  dynamics  (CFD)  technique  is  introduced 
for  the  stack  modelling.  More  specifically,  the  com¬ 
mercial  CFD  package  PHOENICS  [6]  is  implemented 
with  formulations  which  describe  the  physical  and  chem¬ 
ical  processes  in  a  cell.  The  stationary  results,  such  as 
the  three-dimensional  distributions  of  temperature, 
pressure,  density  and  gas  composition  across  the  cell 
or  stack,  are  presented  as  examples  of  the  model 
applications.  Further  validation  of  this  model  and  the 
study  of  the  dynamic  performances  with  regard  to  load¬ 
following  mode  will  be  continued. 


2.  Description  of  main  processes 

Fuel  cells  convert  chemical  energy  directly  into  elec¬ 
trical  energy.  The  principle  of  MCFC  operation  is  shown 
in  Fig.  1.  The  process  involves  the  oxidation  of  a  fuel 
(such  as  hydrogen)  and  the  reduction  of  an  oxidant 
(usually  oxygen).  Electrons  pass  from  the  fuel  electrode 
to  the  oxidant  electrode  via  an  external  load,  and  the 
electrical  circuit  is  completed  by  ionized  particles  cross¬ 
ing  an  electrolyte.  The  main  processes,  which  determine 
the  fuel-cell  stack  performance  with  regard  to  electricity 
generation,  are  identified  as  mass  transport,  chemical, 
heat-transfer  processes,  and  voltage-current  relations. 
These  are  discussed  in  detail  as  follows. 


2.1.  Mass-transport  processes 

The  reactant  gases,  fuel  and  gas  mixture  of  air  and 
C02,  are  separately  fed  to  the  anode  and  cathode 
manifolds  in  the  fuel-cell  stack.  At  the  anode  and 
cathode  inlet  manifolds,  the  gases  are  transported  along 
the  stacking  direction.  Further,  in  the  cell  anode  and 
cathode  flow  paths,  the  gases  flow  from  inlet  manifold 


Fig.  1.  Electricity  generation  principle  of  MCFC. 


to  outlet  manifold.  At  the  anode  and  cathode  outlet 
manifolds,  the  gases  exit  along  the  stacking  direction. 

In  the  fuel  cell  flow  path  at  the  cell  level,  absorbed 
C02  and  02  on  the  cathode  side  react  to  form  C032-, 
which  reacts  at  the  anode  side  with  absorbed  H2  to 
form  H20  and  C02,  with  release  of  electrons. 

2.2.  Chemical  processes 

The  electrochemical  reactions  considered  at  the  fuel¬ 
cell  electrodes  are  the  following: 

(i)  at  the  anode: 

H2(g)  +  C032-  (1)  — ■*  H20(g)  +  COz(g)  +  2e~  (1) 

(ii)  at  the  cathode: 

0.5O2(g)  +  COz(g)  +  2e~  — *  C032  “  (1)  (2) 

The  overall  cell  reaction  is: 

H2(g)  +  0.5O2(g)  +  C02(c)(g)  H20(g)  +  C02(A)(g) 

(3) 

where  the  subscripts  C  and  A  corresponds  to  cathode 
and  anode,  respectively. 

On  the  anode  side,  also  the  water-shift  reaction  occurs, 
increasing  the  H2  concentration: 

CO(g)  +  HzO(g)  *=*  C02(g)  +  H2(g)  (4) 


2.3.  Voltage-current  relation 


The  maximum  electrical  work  obtainable  (IFel)  in  a 
fuel  cell  operating  at  constant  temperature  and  pressure 
is  given  by  the  change  of  the  Gibbs  free  energy  (AG) 
of  the  electrochemical  reaction: 


Wel  =  AG  =  -nFE 


(5) 


The  maximum  electrical  voltage  obtainable  (E)  in  a 
single  fuel  cell,  which  is  often  called  reversible  cell 
potential,  is  given  by  the  Nernst  equation: 


E=E°+  —  In 
2 F 


yuj  oWo 
FH2QVc02(A) 


(6) 


For  the  overall  cell  reaction,  the  cell  potential  in¬ 
creases  with  an  increase  in  the  activity  of  reactants 
and  a  decrease  in  the  activity  of  products.  Changes  in 
temperature  also  influence  the  reversible  cell  potential. 

However,  the  working  voltage  (V)  for  a  fuel  cell  is 
much  lower  than  the  Nernst  potential.  There  are  ir¬ 
reversible  losses  in  a  practical  fuel  cell,  which  primarily 
result  from: 


V=E-  AFdiff(g)  -  AFconc  -  AKkine  -  AKdiff(I)  -  AFcond 

(7) 

where  AFdiff(g)  is  the  gas  diffusion,  AFcond  the  electronic 
conductivity,  AKkine  the  electrode  kinetics,  AFdiff(])  the 
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liquid-phase  diffusion,  and  ALcond  the  finite  ionic  con¬ 
ductivity. 

At  the  stack  modelling  level,  there  are  additional 
sources  of  loss  (e.g.  inhomogeneous  gas  flow  distribution 
among  the  cells).  In  the  present  fuel  cell  model,  the 
overpotential  phenomena  are  simply  considered  as  a 
voltage-current  relation. 

2.4.  Heat-transfer  processes 

The  heat  generation  in  a  fuel  cell  is  a  measure  of 
its  inefficiency,  but  it  is  also  very  important  for  ap¬ 
plications  such  as  cogeneration.  It  is  helpful  to  distin¬ 
guish  two  types  of  heat  sources.  One  is  the  heat 
associated  with  the  reversible  cell  reaction.  The  other 
is  the  irreversible  heat,  caused  by  the  resistance  to 
current  passage.  In  conventional  stacks,  operating  on 
externally  reformed  natural  gas,  significant  removal  of 
heat  is  necessary  to  maintain  the  stack  temperature  at 
about  650  °C. 

There  are  conductive,  convective  and  radiative  heat- 
transfer  phenomena  between  the  fuel  gas  at  the  anode 
side  compartment,  electrodes  with  electrolyte  matrix 
and  oxidant  gas  at  the  cathode  side  compartment  in 
the  fuel  cell.  Furthermore,  there  is  heat  transfer  due 
to  mass  transfer  between  the  fuel  gas  and  oxidant  gas. 


3.  Related  physical  laws 

The  conservation  laws  of  mass,  energy  and  momentum 
are  the  basis  of  studying  fuel  cells  [7].  Furthermore, 
electrochemical  equations  such  as  Nernst  equation  (Eq. 
(6)),  give  the  electrical  performance  of  fuel  cells  ac¬ 
cording  to  the  energy  conservation. 

The  terms  in  the  following  differential  equations 
denote  influences  on  a  unit-volume  basis.  Thus: 

a jx  a l  a l 

net  efflux  per  unit  volume  =  div  /= - 1 — z  H - 

<k  ty  fe 

(8) 

3.1.  Conservation  of  a  chemical  species 

Let  m,  denote  the  mass  fraction  of  a  chemical  species. 
In  the  presence  of  a  velocity  field  u,  the  conservation 
of  m,  is  expressed  as: 

a 

—  (pm,)  +  di  \(pum,  +/,)  =R,  (9) 

ot 

Here  d(pm,)/dt  denotes  the  rate  of  change  of  the 
mass  of  the  chemical  species  per  unit  volume.  The 
quantity  pum,  is  the  convection  flux  of  the  species,  i.e., 
the  flux  carried  by  the  general  flow  field  pu.  The  symbol 
J,  stands  for  the  diffusion  flux,  which  is  normally  caused 


by  the  gradients  of  m,.  The  quantity  R,  on  the  right- 
hand  side  is  the  rate  of  generation  of  the  chemical 
species  per  unit  volume.  This  generation  is  caused  by 
the  chemical  reaction. 

3.2.  Energy  equation 


The  energy  equation  can  be  written  as: 


where  h  is  the  specific  enthalpy,  k  the  thermal  con¬ 
ductivity,  and  Sh  the  volumetric  rate  of  heat  generation. 


3.3.  Momentum  equation 

The  differential  equation  governing  the  conservation 
of  momentum  in  a  given  direction  for  a  newtonian 
fluid  and  laminar  flow,  which  are  assumed  in  this  paper, 
can  be  written  as: 

+div(puu- g.Vu)=  -  ^  (11) 

where  p.  is  the  viscosity,  and  p  the  pressure. 


4.  Computational  fluid  dynamic  approach 

The  software  package  PHOENICS  provides  a  solver 
for  fluid-flow,  heat-transfer,  chemical-reaction  and  re¬ 
lated  problems  ranging  in  complexity  from  one-dimen¬ 
sional  single  phase  and  steady  state,  to  three-dimen¬ 
sional  multi  phase  and  transient.  PHOENICS  is  used 
here  to  solve  the  three-dimensional  conservation  equa¬ 
tions  (Eqs.  8  to  11),  which  have  been  described  in  the 
former  sections,  in  finite-volume  form.  All  these  con¬ 
servation  equations  can  be  expressed  in  a  single  form: 

|  (pf)  +  div(pF<f>  -  /;grad  <£)=S*  (12) 

where  r*  is  the  diffusion  coefficient,  S#  the  source 
term,  and  stands  for  any  one  of  the  dependent 
quantities:  m„  h,  u  or  /.  When  4>—l,  Eq.  (12)  turns 
to  be  a  continuity  equation. 

To  close  the  equations  for  a  unique  solution,  boundary 
conditions  must  be  provided.  A  universal  formula  can 
be  used  to  describe  all  kinds  of  boundary  conditions 
encountered  in  a  fuel  cell  as  a  source  term  in  Eq.  (12). 

5^=[Q  +  Cmass(Emass-F)](^-<A)  (13) 

where  C*  and  V#  are  the  ‘coefficient’  and  the  ‘value’, 
which  are  prescribed  separately  for  each  variable.  Fmass 
has  the  physical  significance  of  an  external  pressure, 
and  the  coefficient  Cmass,  that  of  the  reciprocal  of  a 
flow  resistance  between  the  internal  grid  node,  where 
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the  pressure  is  p,  and  the  external  pressure  is  Kmass 
[6].  In  this  study,  several  separated  sub-routines  have 
been  developed  for  the  implementation  of  the  chemical 
reaction,  voltage-current  relation  and  physical  prop¬ 
erties  of  gas  and  fuel-cell  hardware.  The  fuel-cell 
modelling  has  been  performed  using  PHOENICS  ver¬ 
sion  1.6.6  on  Convex  3840,  and  the  following  dem¬ 
onstrations  were  conducted. 


5.  Setting-up  cases 

To  evaluate  the  present  fuel-cell  model,  four  cases 
listed  in  Table  1  were  selected  as  examples  for  dem¬ 
onstration.  The  first  three  cell  cases  are  used  to  illustrate 
the  three-dimensional  distributions  of  temperature, 
pressure,  concentration  and  density  across  fuel  cells 

Table  1 

Main  features  of  the  four  cases 


Objectives  of  demonstration 


Ceil  with  cross  flow 


Temperature  profiles  across 
the  cell  at  full  and  part 
electrical  load 


Cell  with  coflow 


Distribution  of  temperature, 
pressure,  concentration  and 
density  along  the  cell 
direction 


with  three  different  manifolds.  By  checking  the  be¬ 
haviour  of  these  distributions,  the  performance  of  the 
model  dealing  with  the  fundamental  processes  can  be 
examined.  Then,  the  stack  case  is  set  up  to  calculate 
temperature  performance  and  mass  distribution  across 
the  stack. 

These  four  cases  are  assumed  and  characterized  by 
the  following  specifications. 

•  Each  cell  has  an  electrochemical  active  area  of 
1  mXl  m  and  a  thickness  of  0.014  m  surrounded  by 
a  frame  of  20  mm  containing  the  internal  gas  manifold. 

•  The  stack  contains  5  cross-flow  cells.  The  inlet  and 
the  outlet  of  the  fuel  and  oxidant  gases  are  located  at 
the  bottom  of  the  stack.  This  means  that  the  fuel  or 
oxidant  gas  first  goes  along  the  stack  direction  within 
the  inlet  manifolds,  then  turns  along  the  cell  flow  path, 
and  finally  leaves  downwards  the  stack  via  the  outlet 
manifolds. 

•  The  fuel  gas  uses  Dutch  natural  gas  (a  mixture 
gas  of  H2,  CH4,  N2,  02,  CO,  C02,  HzO)  which  is 
supplied  to  the  anode  from  the  external  reformer  [8], 
and  the  oxidant  gas  uses  a  mixture  gas  of  air  and  C02 
which  is  supplied  to  the  cathode. 

•  The  inlet  temperature  of  fuel  and  oxidant  gas  is 
600  °C,  and  the  outlet  is  at  atmosphere  pressure.  Both 
the  fuel  and  oxidant  gases  are  considered  as  ideal. 

•  The  operation  of  cell  at  constant  voltage  is  also 
imposed,  and  average  current  density  is  set. 

•  The  criteria  for  safe  cell  or  stack  operation  here 
only  monitor  its  temperature  and  pressure  performance. 
That  is,  severe  operational  conditions  only  indicate  high 
or  high  gradient  of  temperature  on  cell  or  stack  hard¬ 
ware,  and/or  high  pressure  difference  between  fuel  and 
oxidant  gases. 


6.  Results  and  discussion 


Cell  with  counter  flow 


Distribution  of  temperature 
and  pressure  along  the  cell 
direction 


Stack  with 


flow 


Temperature  distribution 
across  stack  and  mass 
distribution  along  the  stack 
direction 


The  results  of  the  four  cases  will  be  discussed  in 
the  following  sections.  The  height  of  cell  or  stack  has 
been  enlarged  by  a  scale  of  20:1. 

6.1.  Distributions  of  main  parameters  across  fuel  cells 

6.1.1.  Temperature  distribution  with  cross-flow 

In  this  Section  we  will  examine  the  temperature 
distributions  in  the  fuel  cell,  which  can  effectively 
illustrate  the  basic  features  of  the  operational  per¬ 
formance.  The  following  three  factors  have  a  contri¬ 
bution  to  the  temperature  distribution. 

(i)  reversible  heat  generated  from  chemical  reactions; 

(ii)  irreversible  heat  generated  from  electrical  resis¬ 
tance,  and 

(iii)  heat  transfer  to  the  surroundings. 

The  temperature  profile  of  the  cell  with  cross-flow 
operating  at  an  average  current  of  1500  and  750 
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Fig.  2.  Temperature  profile  of  a  cross-flow  cell  at  an  average  current  density  1500  A/m2  (a)  temperature  at  a  transverse  section  of  cell 
hardware,  and  (b)  isothermal  surface  of  temperature  700  and  630  °C. 


(a)  (b) 

Fig.  3.  Temperature  profile  of  a  cross-flow  cell  at  an  average  current  density  of  750  A/m2  (a)  temperature  distribution  at  a  transverse  section 
of  cell  hardware,  and  (b)  isothermal  surface  of  temperature  650  and  610  °C. 


A/m2  are  shown  in  Figs.  2  and  3.  The  main  features 
may  be  identified  as  follows. 

•  The  point  of  maximum  temperature  is  at  710  °C 
and  is  approximately  located  at  the  oxidant  gas  outlet 
and  near  the  fuel  gas  inlet.  This  may  be  analysed 
qualitatively  using  the  three  factors  given  above.  The 
first  reversible  heat  from  the  chemical  reactions  at  this 
location  is  high,  since  the  fuel  gas  is  relatively  con¬ 
centrated.  The  fuel  gas  utilization  is  much  higher  than 
that  of  the  oxidant,  which  implies  that  the  chemical 
reaction  rate  is  more  dependent  on  fuel  gas  composition 
than  that  of  oxidant.  The  second  factor,  irreversible 
heat  from  the  electrical  resistance,  is  determined  by 
the  current  and  electrical  resistance  according  to  Ohm’s 


law.  The  third  factor  results  in  a  higher  temperature, 
since  the  oxidant  gas  at  this  location  has  been  heated 
as  it  passes  the  cell. 

•  As  expected,  the  maximum  temperature  and  the 
temperature  gradient  of  cells  operating  at  a  current 
density  of  1500  A/m2  is  much  higher  than  at  750  AJ 
m2.  Hence,  a  more  severe  temperature  distribution 
appears  at  high  current  density. 

6.1.2.  Operational  characteristics  with  co-flow 

The  distributions  of  temperature,  pressure,  density 
and  gas  composition  for  co-flow  cells  operating  at  an 
average  current  density  of  1500  A/m2  are  discussed 
below.  The  corresponding  results  are  shown  in  Fig.  4. 
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Fig.  4.  Performance  of  a  co-flow  cell  at  an  average  current  density  of  1500  A/m5:  (a)  temperature:  (b)  pressure  (Pa);  (c)  H2  and  C02 
concentration  (mol/mol),  and  (d)  gas  density  (kg/m3). 


The  main  characteristics  can  be  summarized  as 
follows. 

•  The  temperatures  of  gases  and  cell  hardware 
increase  in  the  direction  of  gas  flow,  and  the  increment 
is  larger  near  the  gas  inlet.  This  is  the  result  of  the 
heat  generated  from  the  chemical  reactions  and  from 
the  electrical  resistance.  In  this  case,  the  point  of 
maximum  temperature  is  at  680  °C,  and  is  located 
around  the  anode  gas  outlet.  However,  the  highest 
temperature  gradient  is  near  the  inlet  between  the 
cathode  gas  and  electrode. 

•  The  pressure  of  fuel  and  oxidant  gas  decreases  in 
the  direction  of  the  gas  flow  as  expected.  The  highest 
pressure  difference  is  4  Pa,  which  occurs  between  the 
inlet  of  the  anode  and  cathode  gases.  It  is  not  large 
because  the  point  of  maximum  pressure  at  the  anode 
corresponds  to  the  point  of  maximum  pressure  at  the 
cathode. 

•  As  expected,  H2  at  the  anode  and  COz  at  the 
cathode  decrease  in  the  direction  of  gas  flow,  and  the 


decrease  being  more  rapid  near  the  inlet  of  the  gas. 
This  implies  a  greater  rate  of  chemical  reaction  in  the 
incoming  gas  than  in  the  depleted  gas. 

•  The  density  of  the  anode  gas  significantly  increases 
and  the  density  of  the  cathode  gas  decreases,  in  the 
direction  of  gas  flow.  In  this  case,  the  density  of  the 
anode  gas  at  the  inlet  is  lower  than  0.17  kg/m3  and  at 
the  outlet  is  greater  than  0.32  kg/m3,  while  the  density 
of  the  cathode  gas  at  the  inlet  is  greater  than  0.42  kg/ 
m3  and  at  the  outlet  it  is  lower  than  0.38  kg/m3.  The 
change  of  gas  density  results  from  three  factors:  com¬ 
position,  temperature,  and  pressure  changes.  A  greater 
percentage  of  the  gas  has  heavier  molecules,  and  is  at 
a  lower  temperature  and  a  higher  pressure,  resulting 
in  a  higher  density.  The  considerable  change  in  gas 
density  shows  the  usual  assumptions  of  constant  gas 
density  in  the  most  of  present  modelling  at  the  cell  or 
stack  level  is  unrealistic. 
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Fig.  5.  Temperature  profile  of  a  counter-flow  cell  at  an  average  current  density  of  1500  A/m2  (a)  temperature,  and  (b)  pressure  (Pa). 


Table  2 

Performance  of  cells  with  three  different  manifolds 


Maximum  temperature 

Largest  pressure  difference 

and  its  location 

and  its  location 

Cross  flow 

710  °C,  at  the  oxidant 

6  Pa,  near  the  oxidant 

gas  outlet  and  near 

inlet  and  the  fuel  gas 

the  fuel  gas  inlet 

outlet  (not  shown) 

Coflow 

680  °C,  around  the 

4  Pa,  between  the  fuel 

fuel  gas  outlet 

and  oxidant  gas  inlets 

Counter  flow 

700  °C,  near  the  fuel 

7  Pa,  between  the  fuel 

gas  inlet 

outlet  and  oxidant  gas 
inlet 

6.1.3.  Temperature  and  pressure  performance  of  a  fuel 
cell  with  counter  flow 

The  temperature  and  pressure  distributions  for  a 
counter  flow  cell  operating  at  an  average  current  density 
of  1500  A/m2  are  illustrated  in  Fig.  5.  The  main  char¬ 
acteristics  may  be  identified  as  follows. 

•  The  point  of  maximum  cell  temperature  is  700  °C, 
located  near  the  inlet  of  the  anode  gas.  The  temperature 
gradient  is  more  uniform  than  that  of  the  co-flow  cell. 

•  The  highest  pressure  difference  is  7  Pa,  which 
occurs  between  the  fuel  gas  outlet  and  the  oxidant  gas 
inlet. 

6.2.  Comparison  of  fuel  cells  with  different  manifolds 

A  comparison  of  fuel  cell  types  may  be  carried  out 
in  various  ways,  depending  on  use  requirements.  Here, 
one  example  is  given  to  illustrate  the  effectiveness  of 
the  model  in  this  application. 

To  determine  the  features  for  safe  and  efficient 
operation  of  fuel  cells,  the  maximum  temperature  of 


the  fuel-cell  hardware  and  the  largest  pressure  difference 
between  fuel  and  oxidant  gases  are  often  considered 
to  be  the  two  most  important  criteria.  Accordingly,  the 
characteristics  of  temperature  and  pressure  difference 
for  fuel  cells  with  different  manifolds  described  in  the 
previous  section  are  summarized  in  Table  2.  The  results 
indicate  that  the  highest  temperature  appears  in  the 
cross-flow  case,  and  the  largest  pressure  difference 
appears  in  the  counter-flow  case. 


6.3.  Performance  across  the  fuel-cell  stack 


With  the  model  described,  the  three-dimensional 
distribution  of  operating  parameters  across  a  stack  may 
be  obtained.  Below,  temperature  profiles  and  mass 
distributions  are  calculated  to  demonstrate  the  use¬ 
fulness  of  the  model  as  a  study  tool. 

The  temperature  distributions  in  a  transverse  section 
of  the  oxidant  gas  outlet  and  the  fuel  gas  outlet  is 
shown  in  Fig.  6,  at  an  average  current  density  of  1500 
A/m2.  The  result  indicates  that  the  highest  temperature 
is  located  at  the  top  cell. 

Mass  flow  along  the  vertical  stack  direction  is  one 
of  the  key  parameters  in  the  stack  model,  since  it  is 
closely  related  to  inhomogeneous  performance  among 
different  cells,  which  may  influence  the  safe  and  efficient 
operation  of  a  stack.  Fig.  7  shows  the  fuel  gas  flow 
distribution  along  the  stack  direction  with  a  cross-flow 
configuration  at  an  average  current  density  of  0,  750 
and  1500  A/m2. 

When  no  current  is  drawn  from  the  fuel-cell  stack, 
the  mass-transport  processes  can  be  modelled  as  a 
network  of  hydraulic  resistances  [9].  Thus,  higher  cells 
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Fig.  6.  Temperature  profile  of  a  cross-flow  stack  of  a  current  density 
1500  A/m2. 


Fig.  7.  Fuel  gas  distribution  along  the  fuel  stack  direction;  current 
density  =  0,  750  and  1500  A/m2. 


(No.  5  being  the  top  cell)  in  Fig.  7  have  less  fuel  gas. 
However,  when  the  current  is  drawn  from  the  fuel-cell 
stack,  the  higher  cells  have  more  fuel  gas.  This  may 
be  interpreted  due  to  the  significant  change  in  fuel- 
gas  density  (see  Fig.  4(d)).  However,  this  requires  further 
confirmation.  From  the  above  results,  it  is  clear  that 


the  model  presented  here  is  capable  of  calculating  the 
mass  distribution  along  the  stack  direction. 

7.  Conclusions 

A  numerical  model  using  the  CFD  technique  for 
analysing  MCFC  performance  is  presented.  The  model 
considers  the  major  factors  which  concern  safe  and 
efficient  electricity  generation.  It  has  been  applied  to 
study  the  three-dimensional  distributions  of  tempera¬ 
ture,  pressure,  gas  concentration,  and  density  across 
fuel  cells  with  different  manifolds.  The  results  appear 
to  be  reasonable.  Moreover,  the  model  has  demonstrated 
features  which  may  be  used  to  optimize  the  design  and 
to  analyse  the  operation  of  fuel  cell  stacks. 

The  conclusions  only  apply  to  internal  manifolding. 
External  manifolding  may  be  different,  e.g.  in  the  cross- 
flows  case. 
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